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Abstract. Usually tunneling is established after imposing some matching 
conditions on the (time-independent) wave function and its first derivative at 
the boundaries of a barrier. Here an alternative scheme is proposed to determine 
tunneling and estimate transmission probabilities in time-dependent problems, 
which takes advantage of the trajectory picture provided by Bohmian mechanics. 
From this theory a general functional expression for the transmission probability 
in terms of the system initial state can be reached. This expression is used here 
to analyze tunneling properties and estimate transmissions in the case of initial 
Gaussian wave packets colliding with ramp-like barriers. 



PACS numbers: 03.75.Xp, 07.79.Cz, 03.65.Xp, 03.65.Ta, 82.20.Xr 
1. Introduction 

Quantum tunneling, the possibility for a system to pass from one state A to another 
state B through an energetic barrier, can be considered as one of the characterizing 
phenomena of quantum mechanics — although it is more general, appearing whenever 
a system is described by a wave equation and there is coupling among evanescent 
waves [1], as it happens in optics [2], for example. In 1928, shortly after the formulation 
of Schrodinger's equation, tunneling was proposed as the physical mechanism that 
explained both field electron emission [3] and alpha decay [4,5]. For about 30 years, 
no satisfactory explanation was possible for these effects, well-known since the end of 
the XlXth century [6,7]. Nowadays tunneling not only appears in tunnel microscopy 
or nuclear physics, but also in other fields, such as semiconductors or catalytic and 
enzymatic reactions, with important direct applications. There is so much written 
about tunneling that it seems there is few space left for new conceptual ideas about 
this phenomenon or mechanism. Nonetheless, it continues stirring our interest. 

In the study and analysis of tunneling problems, many different techniques 
(classical, semiclassical and quantum-mechanical) have been considered in the 
literature [8,9]. For example, for low tunneling (i.e., near the top of the barrier), 
classical and semiclassical approaches, such as the WKB approximation, have been 
shown to be very appropriate both computationally and interpretively. Actually, 
after some refinements deep tunneling (i.e., well below the top of the barrier) can 
also be studied [10]. However, for some other cases, this methodology fails and the 
problem has to be solved exactly quantum- mechanically. In these cases, though, 
interpretations are often based on classical and semiclassical arguments. A way 
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to avoid these thought schemes and provide ahernative arguments fuUy based on 
quantum-mechanical grounds consists of considering Bohmian mechanics [11,12]. The 
first work in this direction was developed by Hirschfelder et al. [13] in 1974, where 
tunneling across a two-dimensional square barrier was studied in terms of stationary 
quantum trajectories or quantum streamlines. As it was shown, under tunneling 
conditions quantum trajectories present analogous behaviors to those observed in 
optics in situations such as the frustrated total reflection or the Goos-Hanchen shift. 
Later on, in 1982, Dewdney and Hiley [14] analyzed the problem of time-dependent 
scattering off square barriers and wells, formerly considered by Goldberg et al. [15] with 
wave packets. More recently, Lopreore and Wyatt started [16, 17] the development of 
the so-called quantum trajectory methods [18] by studying tunneling through barriers 
in one and two dimensions. Within this methodological schemes, quantum-mechanical 
problems are tackled by treating systems as a quantum fluid [19] and then solving the 
corresponding equations of motion as in classical fluid dynamics. The wave function 
or any related property is then synthesized from trajectory calculations, thus skipping 
a direct use of the time-dependent Schrodinger equation. 

An appealing feature to consider Bohmian mechanics in the study of tunneling 
dynamics comes from the fact that initial conditions leading to tunneling can be 
unambiguously determined [16, 17]. This information can be obtained in a rather 
practical way. For example, it allows us to determine how many trajectories pass 
the barrier even if at later times they eventually recross it back and forth again 
(e.g., problems of chemical reaction dynamics [20]). Moved by this fact, here 
we present an alternative scheme to determine tunneling conditions and estimate 
transmission probabilities based on Bohmian mechanics. Accordingly, a connection 
can be established between the transmittance and the system initial state. Thus, first 
we have considered the time-dependent collision of a Gaussian wave packet with a 
ramp-like potential, which can be solved analytically. Then, from this model, tunneling 
properties have been inferred for ramp-like barriertU, such as those we may find in 
problems involving electric fields (e.g., electron ionization processes). In particular, a 
connection between the event-to-event fraction of particles that overcome the barrier 
and the features characterizing both the initial wave packet (width and momentum) 
and the barrier (slope) has been determined. Thus, although quantum trajectories are 
not experimentally observabl^, the information extracted from them could be used in 
a practical fashion in the design or characterization of quantum control experiments 
involving tunneling (e.g., ionization processes or chemical reactions). In this regard, 
it is worth mentioning that this idea has also been considered in the literature to 
determine tunneling times [23] or escape rates for confined multiparticle systems [24] . 

This work is organized as follows. In section [21 the scheme and hypothesis 
considered to infer tunneling information from Bohmian trajectories are discussed. In 
section 131 a practical analysis of this information to estimate transmission probabilities 
and its check against exact numerical calculations are presented and discussed. Finally, 
in section|4l the main conclusions extracted from this work as well as its generalization 
and potential interest in more complex problems are summarized. 

X In this work, we have only focused on transmission below the barrier height and therefore important 
problems like above-barrier reflection are not considered. Work related is currently in progress. 
§ Recently photon paths [21] have been experimentally inferred [22] by using weak measurement 
processes. 
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Figure 1. Schematics of tiie scattering of a Gaussian wave pacltet witti a linear 
ramp potential. If this potential is truncated at some position Xcutoff (i-e., 
V(x > aioutoff) = 0), it is shown that tunneling may take place depending basically 
on how the initial wave packet is prepared and the uphill slope of the resulting 
barrier (see section r2.3l l. 



2. Tunneling with quantum trajectories 

Square barriers are often considered to carry out time- independent analysis of 
tunneling [25,26]. However, the collision of a Gaussian wave packet with a ramp- 
like barrier can be more easily associated with realistic tunneling problems, where the 
ramp represents to some extent the gradual potential slope met by the wave packet. 
Thus, apart from containing the physical elements necessary to understand quantum 
processes and phenomena where tunneling is involved, it is also simple enough to be 
analytically handled. In this regard, consider the collision of a Gaussian wave packet 
with a linear ramp-like potential, V{x) — max, as illustrated in figure [TJ At i = 0, 
the wave packet is given by 

where Aq = {2iTal)-^/'^ is the normalization constant, xq and po ^I'e respectively the 
(initial) position and momentum of its centroid (i.e., {x) = xq and {p) = po = mvo, 
with vo > 0), and ctq its initial spreading. The ramp has a positive slope along 
the x-direction (i.e., a > 0, since the mass m is always positive). Thus, after 
colliding with the potential, the wave packet will move backwards with an acceleration 
a = —{dV/dx)/m ~ —a < 0. However, if the ramp is truncated at some point 
X — Xcutoff (vertical dashed line in figure [Ij , tunneling may appear under certain 
conditions. The purpose below is to determine parameters and conditions leading to 
this tunneling. 

2.1. Scattering with a linear ramp potential 

As shown elsewhere [27], the time-evolution of a free Gaussian wave packet like (HJ 
is governed by two dynamical processes: propagation in configuration space and 
spreading. This fact becomes apparent through the two terms that appear when 
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computing the corresponding average energy or energy expectation value, 

E={H) = ^ + ^,^E,^ (2) 

where H — /2m — —{h?/2m)d'^/dx'^ and Ek stands for the free wave packet 
average kinetic energy. Accordingly, a spreading rate or velocity is defined as 
Vs = Vsl'm = h/2maf), with ps being the associated spreading momentum. The 
dynamical evolution of ([T]) or any other kind of wave packet regardless of its initial 
shape (unless it is spreadless [28-30]) is ultimately ruled by the ratio vo/vs. This can 
be easily seen in diffraction problems, where the "shape" of the initial wave function 
governs its eventual evolution [31-33]. 

In tunneling processes, the ratio vq/vs is also very important. The time-evolution 
of the Gaussian wave packet ([Ij colliding with the linear ramp potential of figure [T] 
can be straightforwardly obtained analytically [34, 35] to yield 

\I>(x t) = j^^g~(^~^<!i)'^/^^t''o+iPcl{x-Xci)/h+i{pl/2m-maxo)t/h-i{po—mat/3)at^/h 

_ ^-{x-x^i)^ /4dtcro+ilPdix-pot/2m)-poxo-ma^t^ /6]/h ^g-j 

where At — {27raf)^^^^ and at = (To(1 + iht/2maQ). The wave packet ^ evolves in 
time without changing its Gaussian shape, but only increasing its width as 




at^\at\^a„\ 1+ =ao\l+^ . (4) 



Its propagation follows the classical uniform accelerated motion displayed by its 
centroid, Xd — xq + v^t — at^/2 and pci = Pa ^ mat. As for its average energy, 

E^Ek + V = ■^+max,i + - ^^E,i + E,, (5) 

Zm smtTQ 

where V = maxg- As it is apparent, this expression contains the classical-like term Ed 
accounting for the translation motion and the (quantum-like) spreading energy. Eg. 
Since there is no coupling between both motions, these two terms remain constant with 
time separately. For more complex potential functions, though, this does not hold; the 
"accommodation" of the wave packet to the corresponding boundary conditions leads 
to their eventual nonseparability. Nevertheless, their initial separability is enough to 
infer the long-time dynamics from the earlier stages of the wave packet evolution [36] . 

Another way to infer the wave packet dynamics is through the probability density 
associated with 

p(x,i) = ^i=e-(----)^/2'^?. (6) 



TTcr, 



The centroid of this function evolves along Xci, while its width spreads at a rate 



dat vlt 



dt at 

Accordingly, ([6]) decelerates at a rate —a as it approaches the ramp and accelerates 
later on again at a rate a after bouncing backwards. The turning point (tp) for 
its centroid (the position where Vd = and the direction of motion changes) is 
x^p = X() + vl/2a, reached at tf^ — vo/a. From ([7]), we also find that at short 
times, the spreading increases uniformly, as ~ (w^ /o'o)t, which indicates an acceleration 
in the wave packet expansion or "boost phase" [16]. Hence, a spreading or boost 
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acceleration, ag = vI/cfq, can be defined, which is different from the dynamical one, 
a. In the short-time regime, the wave packet thus spreads as crj ~ cto + 
i.e., undergoing a classical-like uniformly accelerated expansion. On the contrary, in 
the long-time regime, the spreading rate is constant and equal to w^, which means a 
uniform spreading of the wave packet, at ~ Vgt. 

The Bohmian trajectories associated with (jS]) can be readily obtained after 
integrating the (quantum) equation of motion [11,12] 

1 dS hH , 

X^——=V^\+—r-^^[X-Xc\), (8) 

m ox Am^a^Ot 
where S(x,t) = [hjli) ln[vl'(a;,t)/**(a;,i)] is the real-valued phase of This yields 

=Xci + — 5o, (9) 

with (5o = 2^(0) — xq being the distance between the initial condition of a quantum 
trajectory and the centroid initial position. According to (jH]), the relative distance 
between any two quantum trajectories x\ and xi increases with time as 



X2{t) - Xi{t) 




(10) 



a;2(0)-a;i(0) 

That is, in agreement with ([7]), initially (i.e., short time-scales) the distance between 
trajectories increases quadratically with time f ([TO|) goes like ^ agt^ /2ao), displaying 
later on a speed-up provoked by the boost acceleration. Nonetheless, at even lager 
times, the trajectories eventually undergo a slowed-down linear increase with time 
(-- Vst/ao). 



2.2. Calculation of transmission probabilities with quantum trajectories 

In standard treatments of tunneling [25,26], transmissions are usually determined from 
matching conditions that the wave function and its first derivative have to satisfy at the 
barrier edges. However, in a general, non-analytical tunneling problem, transmission 
probabilities are obtained by computing the quantity 

/•oo 

Too = lim T{t) = lim / p{x,t)dx (11) 

>-oo >oo / 

where T{t) is the restricted probability [37] in the region behind Xcutoff- From a 
Bohmian viewpoint, (jlip is statistically computed by considering only those quantum 
trajectories that reach the region behind the barrier (x > Xcutoff) or transmission 
region (i.e., as in a standard classical Monte-Carlo sampling). The initial conditions are 
randomly distributed according to po{x) = |^'o(a:)|^ and therefore, at any subsequent 
time, the transmitted trajectories will also be randomly distributed according to the 
transmitted probability density pT{x,t) = p{x > Xcutoff,0- 

Consider x"^ denotes the asymptotic position of the first trajectory that 
penetrates into the transmission region and a;J° the position of the last one. These 
trajectories can then be identified with the conditions that define the beginning and 
the end of the transmitted part of the wave function at i — > cxi. This allows us to 
reexpress (fTTj) as 

Too = / ' PTix°°)dx^ = lim / '^ \Tixit))dxit), (12) 
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where x°° denotes the position of a quantum trajectory at t — > oo and confined between 
xf and (i.e., inside the transmission region), p-j-{x°°) is the probability density 
evaluated on x°° and dx°° is the distance (at t — )■ oo) between x°° and a the closest 
neighbor. The same holds for the expression after the second equality, which stresses 
the time-dependence of the quantum trajectories and the evaluation of the probability 
density along them (this is why we have considered explicitly p{x{t)) instead of p{x, t), 
which simply expresses the probability density evaluated at a point a; at a time t). 
This constitutes a very important step. According to Bohmian mechanics, one can 
follow the trajectories backwards in time until reaching their initial conditions, i.e., 
— >■ a;i(0) and x^ — >■ 2:2(0) as t goes from 00 to 0. Therefore, (fT2]) can be recast as 
an integral over initial conditions or, in other words, the initial state of the system, as 
/• 2:2(0) 

Too - / p{x{0))dx{0). (13) 
"'2:1(0) 

From this general result, we find that the knowledge of both the initial state and 
the initial conditions a;i(0) and a;2(0), one can already determine the transmission 
probability (or get an estimate of it) performing the whole dynamical calculation. 

In order to illustrate the previous assertion, consider the Gaussian wave packet 
above. Substituting ^ into 

p{x{t))^^p{xm, (14) 

and, after assuming small differentials in ([6]), 

dx{t) = — dx{Q). (15) 

Then, substituting now (fT4|) and (fTS)) into (|12l) . we obtain (fT3)) . This is a rather simple 
way to analytically proof the validity of this result for the case we are interested 
in here, namely a Gaussian wave packet (it will be further worked out later on). 
However, let us stress that ()13|) it is a general result, which follows from combining the 
quantum continuity equation and the hydrodynamical picture provided by Bohmian 
mechanics, and has also a strong connection with the so-called Born rule of quantum 
mechanics [38] (to some extent, (|13l) could be considered as a consequence of including 
time-dependence in Born's rule). Furthermore, despite this physical argumentation, 
though, a simple analytical proof can also be found which justifies Consider the 
Jacobian 

from which (jl4p and (1151) can be expressed in a more general way as 

p{x{t))^\J\p{x{Q)) and dx{Q) ^ \J\dx{t), (17) 

respectively. These expressions stress the fact that, as in classical mechanics, in 
Bohmian mechanics there is also a causal connection (mapping) between two points 
x{G) and x{t) in configuration space (in classical mechanics this connection is in phase 
space). Therefore, according to a Liouvillian (conservative) viewpoint, a swarm of 
initial conditions enclosed within a certain region C of the configuration space will end 
up in another bound region C of this configuration space at a subsequent time, without 

I Numerical calculations with general wave packets and barriers have also been carried out to test 
it, although they are not included in this work. 
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any loss ol trajectories during the evolution. Hence, if the integrated probabihty in 
C (i.e., the total number of trajectories within it) at t = ti is Pc{ti) and, at some 
subsequent time ^2 the integrated probabihty in C is Pc'it^), then -Pc'(^2) = Pc{ti) 
if the boundary C results from the causal evolution of C, i.e., C — > C when ti — > t2- 
Similarly, here we have that (|lip and (|13l) are equivalent regardless of the probability 
density considered, i.e., 

/ prix^)dx°° = / prix{0))dx{0), (18) 

Jxf Jxi{0) 

provided xi{0) — ?► xf and a;2(0) — >■ xf as t : oo. 

The result lfT3|) actually holds for any restricted probability and not only for 
tunneling transmissions. For example, in problems such as atom-surface scattering [39] 
or slit diffraction [40], quantum trajectories allow us to determine unambiguously in 
the initial probability distribution the boundaries that separate the contributions that 
lead to each specific feature of the corresponding diffraction pattern. More specifically, 
consider the so-called peak area intensity, which is the integral of the probability 
density covered by a given diffraction peak. In the examples mentioned before, the 
peak area intensity can be obtained directly from the initial wave function once the 
two quantum trajectories that delimit a given diffraction peak are known, just by 
computing the number of (quantum trajectory) initial conditions lying between the 
initial conditions corresponding to those delimiting trajectories. This is exactly the 
same problem we face here with respect to tunneling, although instead of dealing with 
diffraction peak intensities, we have transmission probabilities. 



2.3. Tunneling through a linear ramp barrier 

In order to estimate transmission probabilities from ()13|) in the problem we are 
dealing with here, first the initial conditions for the boundary trajectories have to 
be determined. In principle, it is reasonable to assume that 0:2(0) is such that 
Po(a^2(0)) ~ 0, with po(x{0) > 2:2(0)) — 0. Assuming also that po must vanish at 
distances much closer than a:cutoff, we can take 0:2(0) to infinity for practical purposes. 
On the other hand, a:i(0) will correspond to some a:™'" that constitutes the onset 
of transmission, i.e., any trajectory such that x{0) < a:™'" will bounce backwards, 
thus not passing the barrier. Below, some criteria to estimate a:™'" will be given 
based on quantum trajectory considerations. For now, these assumptions allow us to 
analytically approximate (jl3p by 

where erfc denotes the complementary error function [41]. This is a general result 
provided the initial state is described by a Gaussian wave packet. 

Transmission probabilities can be thus determined from (jl9p if a:™'" is known or, 
alternatively, from ([T3|) by means of a random, Monte-Carlo-like sampling of initial 
conditions, a:i(0), distributed according to po. Proceeding in this way, (|13p becomes 

1 ^ 

Too « J™^;^E'5K°(0) e {xr\^T)), (20) 

i—l 

where (5 = 1 if the initial condition is contained within the considered spatial range 
and otherwise — the sampling according to the initial probability density is denoted 
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by the superscript 'po'- In the hmit of large A/", PU)) approaches p^ . as aheady seen 
elsewhere for chemical reactions and atom-surface scattering processes [20,37]. 

Below, we provide some considerations based on the ramp barrier from figure [T] 
to determine x™'". The cutoff is chosen at a distance such that po{xcMtos) ~ in 
order to ensure there is no probability behind the barrier; below, we will assume 
a^cutoff = + Nao, without loss of generality, where TV is related to some probability 
onset. Thus, although ([9]) itself is not a solution of the corresponding tunneling 
problem, it is very useful in the analysis of tunneling through such barriers in terms 
of quantum trajectories. Depending on the initial translational motion two cases can 
be analyzed. 

2.3.1. Tunneling dynamics without initial translational motion. When wg = 0, 
the wave packet initially rests on the classical turning point and therefore, at any 
subsequent time, it will slide down the ramp. Appearance of tunneling will then be 
ruled by which term is dominant in ([9|), the dynamical down- hill accelerated motion 
induced by the potential or the wave packet spreading. The former leads the wave 
packet to move far away from the barrier, while the latter provokes an anti-downhill 
motion by 'pushing' part of the wave packet opposite to its translation direction. 

For the ramp potential, the centroid classical trajectory will evolve backwards 
(toward x < 0) since the very beginning and, in the long-time regime, all quantum 
trajectories will also move backwards since their dynamics is ruled by a, as seen in 
Moreover, due to the Bohmian non-crossing rule [27], apparent from pO)) . the distance 
among them will increase linearly with time. Thus, the feasibility of tunneling will rely 
on the dynamical behavior undergone by the ensemble of trajectories starting from a;o 
onward (with 5^ > 0) for their path can be eventually closer to Xcutoff and therefore 
it is likely they tunnel if the boost is strong enough. In brief, a quantum trajectory 
will undergo tunneling if its turning point satisfies the condition Xtp > a^cutoff • 

The turning point of a quantum trajectory is given by a vanishing of its associated 
velocity ([8]), i.e., when the translation and spreading terms in this relation cancel 
each other, making the corresponding (quantum) trajectory to bend over and evolve 
backwards. Equating the left-hand side of ^ to zero, substituting ([9|) and wq = into 
its right-hand side, and then solving for t the resulting expression eventually yields 

This turning time is given as a function of three parameters: (1) the slope of the 
potential or, in other words, the dynamical acceleration that it imprints on (quantum) 
particles; (2) the initial spreading of the wave packet, which is a measure of its size at 
the "preparation" instant and, therefore, can be somehow controlled experimentally; 
and (3) the initial distance between the corresponding particle and the center of the 
wave packet, this distance being inconceivable from the viewpoint of standard quantum 
mechanics. 

Given a and croi the question is then whether a particle can cross to the other 
side of the barrier or not. This information can be extracted from as follows. It is 
apparent from (j21|l that, in order to observe a "true" turning point (i.e., the particle 
moves forward for a while before getting backwards), the condition 

^0 > (^^) do = <5S (22) 
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must be fulfilled, where we call (5g the critical distance — equivalently, a;(0) > xd^) = 
Xq + 5^ = xq + {a/as)cr,^. From this position onwards, all trajectories will display a 
turning point different from the origin (see section [3]) , while for any trajectory with 
^0 ^ ^0 its own initial position a;(0) is also its turning point. The trajectory for which 
5q = 5q plays the role of a separatrix or boundary between two different dynamical 
behaviors. Note that Sq grows very rapidly with ao (~ CTq), which means that the 
larger the initial spreading, the further away the position of the turning point. In 
other words, turning points will appear close to xo in cases of fast-boosting initially- 
prepared states, while they will appear far away (even in regions where p sa 0) for 
slow-boosting states. 

The position of the turning point for any trajectory with initial condition 
satisfying is determined by substituting (pij) into which yields 

xtp = xo + -[ — \ +7;{ — ] — = xo + ^— + 7; -■ 23) 

Note that for 5q — 5q we obtain xtp = a;c(0), which is the initial position for 
the boundary trajectory. In order to estimate now which initial positions will 
lead more likely to tunneling, we assume they must be close enough to Xcutoff- 
Thus, equating (|23|) to Xcutoff and the solving the resulting equation for (5o, we 
find (5™'" = (2Q;cro/as)^^^[a;cutoff — xo — (croQ!/2as)]^/^, which means that at least 
all trajectories started from 

mill 2acro j <7q ol 

Xq =Xq + \\ \\Xcv,toS- Xq- — (24) 

V tts V 2 tts 

onward are going to display tunneling through the barrier. Like Xc(0), a;™™ also marks 
the starting point of another type of boundary trajectory: a trajectory that sets up 
a difference between those with turning points going beyond Xcutoff and those that 
cannot reach this point. 



2.3.2. Tunneling dynamics with initial translational motion. From 0, when Vq ^ 
the wave packet moves towards the ramp and then bounces backwards (its centroid 
turns at t^^ = Vo/a, when it is at x^^ = xq + Uo/2a). Again, only those trajectories 
with (5o > will have a chance to cross. Nonetheless, now the initial translational 
velocity introduces a new control parameter into the barrier passage process, as seen 
if we substitute ([9]) into ([8]), with vq ^ 0, and then equate to zero the left-hand side, 
which renders 

= vn-at+ t. (25) 

Unlike the case discussed in Section 12.3.11 here vq ^ Q constitutes a second element 
that can enhance tunneling. Equation (|25|) can be solved for t by using different 
techniques employed to solve quartic [42] and cubic [43] equations and different cases 
can be discussed (depending on the value of the several parameters involved in it). 

First, consider ctq is relatively large and therefore there is a relatively slow boosting 
phase. More specifically, the wave packet spreading is assumed to be negligible 
approximately up to t ~ tjjp, which implies a relatively small boost acceleration 
compared with a. Under this condition, a first-order expansion in time of (j25l) renders 

w Wo - at + ^ t, (26) 
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^tP--(l + -^)=itt + ^^o (28) 



from which we find 

^.P = ^^. (27) 
a — Us — 

Of course, this expression is only valid if at ~ ctq or vq is relatively large (in the limit 
Vq — >■ 0, it does not approach the slow-boosting limit of (|21l) '). Equation (l?7)) can be 
further simplified to 

a \ a (jQ ^ 

if the second term in the denominator of ([27]) is very small in comparison with a 
(e.g., for trajectories very near the center of the wave packet, for which can also be 
small). Under this assumption, according to (j28p . the turning time increases linearly 
with 5q. Furthermore, since crj w (Tq, it can be shown that quantum trajectories will 
be essentially parallel to the centroidal one, 

Xtp K x{Q) + [xf^ - xq) ^ 5o + x%. (29) 

Therefore, for large ctq (slow boosting), we have i5™'" = Xcutoff — a^o — VQ/2a and 
therefore 

x^" = a^cutoff - vl/2a. (30) 

That is, regardless of the value of ctq, there is always an onset of tunneling provided 
the translational velocity is large enough that the trajectory reaches the ramp cutoff. 

The second case is that of fast boosting, i.e., when the initial boosting phase 
occurs much earlier than the trajectory reaches its turning point. Assuming at ~ Vst 
and then substituting this into ([5]) leads to 

« Weff - at, (31) 
where Ves = vq + VsSo/ao is an effective constant velocity encompassing both the 



translational and the spreading velocities, and proportional to Sq. After pip , a 
quantum trajectory reaches its turning point at 

t,^ = ^-f^tll + t's^, (32) 
with t'g^ — (vs/aao)So. As for the turning-point position, we find 

xtp = xo + ^, (33) 

which looks pretty much like the expression for the classical turning point, but with 
Vo substituted by VcS- From ([55]) . we obtain J™'" = (2Q;cro/Q;s)^/^(a;cutoff — xq)^^^ — 
vgao/vs and therefore 



■ 2aao vq 

Xq ^xq + \ va;cutoff - a;o ctq. (34) 

V tts Vs 

The first part in this relation coincides with (|24|) in the limit Xcutoff — a^o S> aoa/2as. 
Hence, in the limit vq — >■ 0, both cases will approach. The interest on ([M)) relies on the 
fact that, unlike (|24p . if the spreading rate is not enough to overcome a;cutoff, there is 
an extra (translational) energy which may help to surmount the barrier, this favoring 
the passage. The starting point of the quantum trajectory can then be further away 
from the barrier cutoff because of the extra velocity. 
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3. Numerical results 

3.1. Parametric study oj the transmission probability 

In figure [5] we show some calculation^ of transmission probabilities, Too^ as a function 
of different parameters in order to illustrate the concepts introduced above. Without 
loss of generality, we have considered vq = 0, which allows us to better understand the 
tunneling dynamics "clean" of possible contributions coming from the kinetic energy 
of the particle. 

Consider first the effect of varying the cutoff distance, a;cutoff- As mentioned 
above, this distance is chosen as a function of the initial width, as Xcutoff =2^0 + Nuq. 
To avoid some arbitrariness in the choice of N and ensure po{x) vanishes far away 
from Xcutoff, we define the ratio 

r(a;eutoff) ^ = 10- (35) 

and call the exponent n the sensitivity parameter^ which provides us with a minimum 
bond for the probability density (relative to its maximum value at the center of the 
initial wave packet). For example, setting n = % means that whenever r(a;(0)) < 
r(xcutoff) ~ 10^®, the probability density will be assumed to be zero (i.e., no 
probability below that value will be detected or computed). Accordingly, we have 
A'^ = VSnln 10 sa 2.15 n^/^. In figure Ufa) we plot Xcutoff (dashed lines) and the 
corresponding curve x™™ (solid lines) as a function of ctq for four different values of 
the sensitivity parameter: n = 4 (green), n = 5 (blue), n = 6 (purple) and n = 7 
(red). As it can be noticed, as n increases, a;™'" also increases, which means to move 
the integration range in (jl3p towards regions with lower values of po{x), this leading to 
smaller transmissions (see figure[2][b)). Nonetheless, the most important changes in the 
integration range happen at high values of ctq , near the onset of no tunneling, this being 
the reason why in figure[21[b) there is only a difference 7^ (n — 4)— 7^(n = 7) ~ 0.0336 
6.72% when this difference is referred to 0.5, the maximum value of Toe) between 
the curves for n — A and n = 7, at ctq ~ 0.138 (see the position of the vertical dashed 
line). In figure [2l[a) we have also plotted (black dashed-dotted line), which does 
not depend on the sensitivity parameter. 

In the central part of figure [2] the effects of the slope of the potential on tunneling 
are displayed for three different values of the dynamical acceleration: a = 5 (green) , 
a = 10 (blue) and a = 20 (red). In all cases, the cutoff corresponds to a sensitivity 
parameter n = 6, with the black solid line of figure[2l[c) denoting the position of Xcutoff 
for each value of ctq. In figure [2jc), for the same value of a, solid lines indicate the 
position of x™™ and dashed ones the position of 5q. As the barrier slope a increases, 
tunneling probability decreases (see figure [DJd)). This decrease is again due to shift of 
the boundary curve 5^ when a increases, which leads to a smaller area between ^cutoff 
and well as to areas where the value of becomes meaningless (see below). 

Taking this fact into account, we can distinguish between two limits of interest. In 
the case of a — > oo, i.e., a vertical, infinite wall, no particle will be able to pass the 
barrier. On the contrary, when a — and we have a flat surface, only half the number 
of particles from the initial ensemble are going to pass the barrier. In this latter case, 
the maximum value of Too is only 0.5 because, as can be easily noted, x^i = and, 

^ In all the calculations presented here, h = m = 1. Moreover, when necessary, quantum trajectories 
have been calculated using standard wave packet propagation techniques [44] to obtain the wave 
function and then integrating l8t . 
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Figure 2. Top: Too as a function of ctq for different values of the sensitivity 
parameter n and, therefore, of the cutoff position x^xxtoS'- n = A (green), n = 5 
(blue), n = 6 (purple) and n = 7 (red). In part (a), scheme to choose the initial 
positions, 2;(0), for the calculation of Too in panel (b): the cutoff position, ^cutoff i 
is indicated by the dashed lines, a:™'" by the solid lines and <5g by the black 
dashed-dotted line. Center: Too as a function of (tq for different values of the 
barrier slope: o = 5 (green), a = 10 (blue) and o = 20 (red). In part (c), scheme 
to choose the initial positions, x(0), for the calculation of Too in panel (d): the 
cutoff position, Xcutoffi is indicated by the black solid line, by the colored 

solid lines and Sq by the colored dashed-dotted line. Bottom: Contributions to 
Too from regions with different sensitivity range (see text for details). In part (e), 
a;™'^ is indicated by the thick solid line, while the straight lines indicate limit of 
regions with sensitivities up to the value given in the corresponding label. In part 
(f), the thick solid line indicates Too, while the colored lines indicate indicate the 
contributions between two consecutive sensitivity values; each color corresponds 
to the upper sensitivity considered from part (e). 
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therefore, half of the particles will move forward (those starting in the front part of 
the wave packet, for which 5q > 0) and another half will move backward (those from 
the rear part of the wave packet, with Sq < 0). 

To understand the different contributions of the region bounded between Xcutoff 
and ccjf™ and therefore their role in the tunneling process, in the bottom part of figure[2] 
Too is plotted as a function of different sensitivity parameters. In figure [2l^e) , we have 
plotted a;™'" (thick black line) for a — 10 and a sensitivity of n = 6 (purple line). 
The positions of Xcutoff for different sensitivities are also plotted. These additional 
lines indicate that if n = 2 (green line), for example, the initial conditions lying 
between this line and x™'" will give rise to probabilities such that r(x(0)) > 10~^. 
The contribution to the transmitted probability of initial conditions lying between 
two consecutive sensitivities n — 1 and n is represented in figure [5Kf); 7^ is denoted 
with thick black line and labels indicate the cutoff chosen, e.g., n — 2 (green line) 
represents 7^ from trajectories starting between the cutoffs with sensitivities n = 1 
and n = 2 (the red and green lines in figure [^Je)). As can be seen, the most important 
contribution (almost 100%) arises from trajectories starting between x™'" and n ~ 1, 
while as n increases the contributions decrease in about one order of magnitude and 
are only relevant as ctq increases. 

3.2. Quantum trajectories and transmission probability 

The results shown above provide us with valuable information concerning tunneling, 
although one can go a step further away yet. In figure |3] different sets of quantum 
trajectories are displayed in terms of ctq- Thus, from top to bottom: CTq^'' = 0.15 (upper 
row), ctq^'' = 0.3 (middle row) and a^^^ — 0.5 (lower row), which correspond (see 
figure [2Ke)) to fast (strong), medium and slow (weak) boosting phases, respectively. 
In all cases, Xq = 0, Vq = and a = 10. Moreover, for computational convenience, we 
have considered a potential function 



where Vq = ma{xQ — 3A^cro) and N — 5.26, which corresponds to n = 6. The left-hand 
side truncation of the ramp avoids the increase of downhill velocity of the wave packet; 
to minimize effects due to the reflections induced by this change of slope [25,26], we 
have assumed a;_ = Xq — SNag. 

In figure [3l[a) we notice that the fast or strong boosting (ctq — ctq^'' ) provokes 
an almost immediate pushing of some trajectories in the opposite direction to the 
evolution of the center of the wave packet (thicker red line). In this case, the expansion 
of the wave packet is such that a.t t — 1 its width has increased more than 20 times (Tq , 
which is more than enough for an important number of trajectories to overcome the 
barrier (see figure [3ljb)). From a standard quantum perspective (see figure [DJc)), the 
very narrow initial wave packet (black) passes to a very wide final one (blue dashed 
line), which splits into a reflected and a transmitted wave packet (blue solid lines) when 
the cutoff is considered (vertical green dashed line). For a weaker boosting (middle 
row), the width of the wave packet at t = 1 is about 5.6cto (see figure[3jd)). Meanwhile, 
the cutoff is at Xcutoff ~ 1.6, which is only 3.5 times smaller than at at t ^ 1. Thus, 
although some trajectories will be able to cross the barrier (see figure |3je)), they start 
in regions where the initial probability density has small values and do not produce 




X < x^ 

X— ^ X -^cutoff 
X ^ a^cutoff 



(36) 
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t t X 



Figure 3. Quantum trajectories illustrating the collision dynamics with a linear, 
ramp-like potential (left-column panels) and the feasibility of tunneling when it is 
truncated (middle-column panels) for: ctq^' = 0.15 (top row), ctq^' = 0.3 (central 

row) and ctq = 0.5 (bottom row). In each trajectory plot, the thick red curve 
indicates the average or expectation value of the position of the wave packet, (x), 
with time. The position of ^'cutoff for the corresponding value of cro is shown by 
the dashed green line (horizontal in the middle-column panels and vertical in the 
right-column ones). In the right-column panels, from top to bottom, the initial 
(black line) and final (blue line) probability densities corresponding to each value 
of fTQ are displayed. In these panels, final probability densities are plotted with 
solid line for the truncated potential and dashed line for the linear, ramp-like one. 



much transmission. As seen in figure [3{f), transmission is almost negligible (about 
500 times smaller than reflection). Finally, in the lower row of figure [31 a case of 
slow or weak boosting where quantum trajectories are out of the conditions leading to 
tunneling is displayed (see figure [2Ke)) — the width of the wave packet here is about 
2.2 times its initial width at t = 1, which is smaller than Xcutoff- This is apparent in 
figure [3fh), where there are not transmitted trajectories, and also in Fig. |31[i), where 
there is no transmitted wave packet. 
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Figure 4. Transmission probability, Too, as a function of the initial width of the 
wave packet estimated with II19I I and II 241 1 (7^"', solid line) and obtained from a 
wave packet simulation {T^^ , circles) for several values of ctq. In the upper inset, 
the same plot but at a logarithmic scale to show the discrepancy between both 
calculations as cro increases. In the lower inset, deviation S (see text for details) 
for the estimated values (black circles) and the corrected ones (red triangles) with 
respect to the values obtained from the wave packet simulation. The values of the 
parameters considered are: a = 10, xo = and vq = 0. 



Finally, in figure 2] we have plotted Too when it is estimated with the previous 
simple model (7^"', solid line) and also from wave packet calculations {T^^, black 
circles) — using quantum trajectories sampled according to po would render exactly 
the same results, as shown elsewhere [39]. As can be noticed, though the estimated 
transmission 7^"' follows very nicely the trend of the correct value, 7^^, there is a 
discrepancy between both magnitudes which increases with ctq. This can be better 
appreciated in the upper inset of the figure, where the Too axis is given in logarithmic 
scale, as well as in the lower inset, where the deviation between both values (measured 
in percentages, %), 



is calculated. In order to understand these deviations, in figures ^e) and [5]Jf) we 
see that, given ctq, those trajectories started between the (red) line for n = 1 and 
a;™'" would be the ones with a larger contribution (the following, between n — I 
and n = 2, contribute with a probability almost two orders of magnitude smaller). 
Thus, note that, beyond ctq sa 0.225 the probabilities are very low (see green curve in 
figure [2](f)) and, therefore, small variations in the position of Xmin may imply larger 
relative errors, as can be seen in the lower inset of figured Now, the source for these 
deviations arises from the fact that, in the model considered here, a;™'" was chosen 
assuming that the last trajectory penetrating into the transmission region coincides 
with Xcutoff- This implies to neglect all those trajectories whose turning points are 
very close to Xcutoff , but without touching it (i.e., with x(0) < a;™'"). Thus, provided 
they are in a neighborhood of Xcutoff, the corresponding trajectories should also be 
taken into account and the a;™'" curve should be correspondingly corrected. If we 
use the trajectory plots to determine the initial position of the boundary trajectory, 
and consider it to estimate again Toa through (fT9)) but with the new x™'" positions, 




(37) 
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effectively, we obtain 7^^, as indicated by the red triangles in the lower panel of 
figure m In order to appreciate the importance of a good characterization of the region 
of initial conditions leading to tunneling, consider, for example, ctq = 0.15. In this 
™, es ^ 0.1774 and 7^^* ~ 0.11795. However, when we localize the initial 
position of the boundary trajectory, x™™' k, 0.1547, and correct accordingly 
the lower value of initial positions, we obtain T^^^ ~ 0.15149, which is very close 
to the value obtained from the simulation, ~ 0.15149. Thus, a displacement 

x'^'''' - X™"' « 0.0227 in an interval xg"'"^ - x™"' w 0.6328 leads to 
a decrease of the true transmission probability of about 22%, while the corrected 
one differs only in about a 0.03% (both value are obtained using (|37l) '). Now, if 

J • •! 1 no -G J rain, est vain, corr r, nono • 

we proceed similarly with ctq = 0.3, we nnd Xq — Xq w 0.0823 m an 

interval Xq^*'°^ — x™™' ~ 0.7024. This leads to a decrease of the true transmission 
probability of about 60%, while the corrected one differs only less than 0.5%. 

4. Conclusions 

Usually in standard quantum mechanics, the transmission or tunneling probability 
is associated with the height and width of a barrier. Analytical expressions for 
this probability can be derived from time- independent or semi-classical calculations 
[25,26] after setting some matching conditions on both the wave function and its first 
derivative at the barrier edges. On the other hand, from a numerical viewpoint, 
transmission probabilities are usually obtained from wave packet calculations, 
computing the amount of transmitted probability accumulated beyond a certain 
boundary [20]. 

A complementary and alternative way to look at tunneling readily arises when 
Bohmian mechanics is considered. In principle, this approach helps us to discern which 
part of the incident wave packet is transmitted through the barrier and how this takes 
place by monitoring the trajectory fiow along time. However, as we have shown here, 
one can go beyond these facts and determine tunneling conditions by studying the 
individual behavior of quantum trajectories, which can be used in a practical way 
to obtain fair estimates of transmission probabilities with information only related 
to the system initial state, according to the general result p^ . To illustrate this 
interesting fact, we have considered the collision of a Gaussian wave packet with a 
ramp-like barrier, for which (|13|) can be recast as in the form given by (|19p . This 
general expression for Gaussian wave packets presents some remarkable properties. 
First, tunneling can be explained essentially in terms of three physical (measurable) 
parameters: the wave packet motion (through vq), the wave packet spreading (through 
at and therefore v^) and the barrier slope (through a). This arises through the 
estimates of x™'", according to expressions like ([M)) . (pO)) or ([34]). Second, once xjf' " 
is set up, (|19p allows us to estimate and compute transmission probabilities. Third, 
conversely, given the transmission probability the inverse procedure can be used to 
determine the dynamical boundaries of the wave packet that lead to tunneling, which 
can be used to extract valuable information about the physical properties associated 
with the wave packet or the barrier (using expressions like (p4|) . (|30l) or (|34|) '). 

Based on the aforementioned properties, in our opinion, the study presented 
here constitutes an important resource at an applied level. Despite that quantum 
trajectories are not experimentally observable, the information they provide can be 
seen as a tool to analyze experimental processes and phenomena where tunneling 
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is involved. In particular, it could be employed to understand and implement 
mechanisms aimed at quantum controlling molecular systems, alternative to (or 
cooperative with) other mechanisms proposed in the literature [45-51], since the 
treatment here described stresses a direct relationship between the experimental effect 
(tunneling) and the initial state. In virtue of this relationship, fairly well summarized 
by (fT3|) or its Gaussian version, (fT9| , one could control the further state of the system 
(i.e., the occurrence of tunneling) by selecting different values of the parameters 
involved in the preparation of the initial state (i.e., the initial conditions of the 
Gaussian wave packet). 
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